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ABSTRACT 

We study the applicability of the idea of internal absorption of 7-rays pro- 
duced through synchrotron radiation of ultrarelativistic protons in highly mag- 
netized blobs to 1ES 0229+200 and 3C 66A, the two TeV blazars which show 
unusually hard intrinsic 7-ray spectra after being corrected for the intergalactic 
absorption. We show that for certain combinations of reasonable model param- 
eters, even with quite modest energy requirements, the scenario allows a self- 
consistent explanation of the non-thermal emission of these objects in the keV, 
GeV, and TeV energy bands. 

Subject headings: Radiation mechanisms: non-thermal - Gamma rays: galaxies - 
X-rays: galaxies - Galaxies: active - BL Lacertae objects: individual: 1ES 0229+200 
- BL Lacertae objects: individual: 3C 66A 
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1. Introduction 

Over the last few years a number of active galactic nuclei (AGN) with redshifts z > 0.1 
have been detected in the very high energy (VHE; E > 100 GeV) regimej The detection 
of VHE 7-rays from such distant objects implies serious constraints on the intensity and 
spectral shape of the extragalactic background light (EBL). Traveling over cosmological 
distances, high energ y 7-rays are significantly absorbed due to effective interactions with 



photons of the EBL (INikishov 



1962 



Gould &: Schreder 



19671 ). The level of attenuation 



depends strongly on the intensity, spectral shape and redshift-depen dence of the EBL. 



Whi 



e robust EBL lower limits can be obtained from galaxy counts ( IMadau Sz Pozzetti 



20001 ). the derivation of the EBL properties based o n direct measurement s is quite diffic ult 



becau se of dominant foregrounds (see for a review 



1998 


). In this regard, the theoretical moc 


eling 


e.g. 


3 rimack et al. 


2008; 


Kneiske & Dole 


2010 


2010 


) is an important aspect of the activity in 



Hauser fc Dwek 



2001 



Hauser et al. 



Franceschini et al. 



2008 



Dominguez et al. 



The mean free path of 7-rays due to interactions with EBL strongly depends on energy. 
Therefore the intergalactic absorption leads not only to attenuation of the absolute fluxes, 
but also to significant changes in the spectral shape of 7-rays. The proper understanding 
of this spectral deformation is crucial for the correct interpretation of the VHE data from 
distant AGN. It is important to note that because of strong Doppler boosting of the 
non-thermal emission (F 7 cx <5 4 ) the 7-ray emission from the brightest blazars can remain 
detectable even after severe intergalactic absorption. 

Thus the mere attenuation of the 7-ray emission is not enough to derive robust 
constraints on the EBL models. The distortion of the initial spectral shape of 7-rays 



1 See \http://tevcat. uchicago. edu/ for an updated list of VHE 7-ray sources 
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contains more information. Since in the effective absorption regime the optical depth 
r > 1, even a slight change of the EBL intensity can lead to a strong change of the 
energy-dependent spectral deformation factor exp[—r(E)]. This allows quite meaningful 
upper limits on the EBL in the relevant energy bands, based on the condition that the 
intrinsic spectrum of 7-rays should have a decent form, e.g. be not much harder than E~ 2 . 
On the other hand, the absorption-corrected VHE spectra of some AGN with z > 0.1 
in some cases appear very hard, even for very a low EBL flux, with a power-law photon 
i ndex T\ nt (dN/dE = Nt)E~ r i M ) quite close to the har dest conventional value of r int = 1.5 



(lAharonian et al 



2006 



2007 



Franceschini et al 



20081 ). In the case of slightly higher fluxes 



of the EBL, the reconstructed spectra would get even harder, with Ti nt < 1.5. 

Although currently there is a general consensus in the community that the EBL 
intensity should be quite close to the robust lower limits derived from galaxy counts, the 
possibility of slightly higher fluxes of the EBL cannot yet be exclude d. In particular, 



using 



Spitzer data and a profile fitting of the faint fringes of galaxies, 



Levenson fc Wright 



(120081 ) claimed a new fiducial value for the contribution of galaxies to the EBL at 3.6 //m 



of 9.0 + ^iiWm- 2 sr\ w! 



rich exceeds 



Franceschini et al. 



jy a factor of ~ 1.6 the flu x of the EBL suggested 



(12008 ). Following 



Levenson fc Wright 



(1200a ] 



Krennrich et al. 



fl2QQ8h 



by 



indicated that for this flux of EBL the initial (absorption corrected) VHE spectra of distant 
blazars 1ES 0229+200, 1ES 1218+30.4 and 1ES 1101-232 (located at redshifts z = 0.1396, 
0.182 and 0.186, respectively) would have a photon index < 1.3. This result would challenge 
the conventional models for VHE production in AGN. 

Generally, the X- and 7-ray non-thermal emission of blazars is interpreted as a sum of 
synchrotron and inverse Compton (IC) components of radiation from relativistic electrons, 
in the framework of the so-called synchrotron self-Compton (SSC) or external Compton 
(EC) scenarios. In the case of radiatively efficient models, i.e. assuming a radiatively cooled 
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particle distribution, the IC spectrum in the Thomson limit is expected to be steeper 
than the power- law distribution with photon index 1.5. This limit does not depend on 
the electron initial (injection) spectrum and can be achieved, for example, in the case of 
a mono-energetic injection. At higher energies, the 7-ray spectrum becomes steeper due 
to the Klein-Nishina effect. We note however, that typically the spectra obtained in the 
frameworks of SSC scenario are steeper, with photon indices ~ 2. Therefore, the spectrum 
with photon index Ti nt = 1.5 is often referred to as the hardest spectrum allowed by 
standard blazar models. However, in the expense of radiation efficiency it is possible to 
produce harder VHE spectra still within the SSC framework, for exam ple assuming a high 



lower-energy cutoff in the electron spectrum (jKatarzyriski et al. 



20061 ). The postulation of 



such a cutoff in the electron spectrum implies very low efficiency of radiative cooling which, 
in turn, increases the requirements for the energy in accelerated electrons and at the same 
time requires very small magnetic fields. Thus, in such scenarios we face a significant (by 



order s of magnitude) deviation from equipartition, W c » Wb (see e.g. 



Tavecchio et al. 



200J). 



Alternatively, 



Aharonian et al. 



( 120081 ) have suggested a scenario for the formation of 
VHE spectra of almost arbitrary hardness by involving additional absorption of VHE 7-rays 
interacting with dense radiation fields in the vicinity of the 7-ray production region. The 
key element in this scenario is the presence of a dense photon field with a narrow energy 
distribution or with a sharp low energy cut-off around > 10 eV. In this case, 7-rays are 
attenuated more effectively at energies ~ 100 GeV than at energies ~ 1 — lOTeV, and 
therefore, for large optical depths (r > 1), the emerging spectrum in the VHE band should 



gradually harden towards higher energies (for detail, see 



Aharonian et al. 



2008|). 



While the absorption of high energy 7-r ays in the inner 



possible, or even unavoidable in some cases ( iMcBreen 



1979 



parts of AGN jets is generally 



Liu fc Bai 



2006 



Reimer 



20071 : 
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Sitarek fc Bednarekl 120081 : iLiu et alJ 12008 Bai et al.ll200a iTavecchio fe Mazinl 120091 ). the 
detailed modeling of this process requires additional assumptions concerning the presence 
of low-frequency radiation fields, the location and size of the 7-ray production region, the 
Doppler factor of the jet, etc. We note that curre ntly there is no observa tional evidence 



Aharonian et al. 



(120081), also in the 



excluding the photon field properties required by 
case of BL Lacs. Remarkably, the internal absorption hypothesis provides an alternative 
explanation for the non-thermal X-r ay emission, namely as synchrotron radiation of 



secondary (pair-produced) electrons (I Aharonian et al. 



20081 ). which suggests a possible 



solution to the pro blem of low acceleration 



emission of blazars ( iCostamante et al. 



20091 ). 



efficiency in leptonic models of high energy 



In the original paper, 



Aharonian et al. 



( 120081 ) presented a general description of the 



scenario with calculations of model SEDs, but the obtained spectra were not compared 
with available data. In the present paper, we discuss the multiwavelength properties of the 
radiation in the internal absorption scenario, and apply the model to the data of two distant 



AGN, namely 1ES 0229+200 



z = .1396) and 3C 66A (estimated at z = 0.444), detected . 



in TeV band (Aharonian et al. 



2007 



Aliu et al. 



2009 



Acciari et al. 



2009 



Reyes et al. 



2009b- 



Here we adopt the proton synchrotron radiation as the source of primary 7-rays, and 
consider the absorption due to 7-7 pair production both in the 7-ray production region and 
in the surroundings. The synchrotron radiation of secondary pairs gives rise to an additional 
lower energy non-thermal component. The latter can be calculated self-consistently and 
depends on the primary 7-ray spectrum, the target photon field and the relativistic motion 
of the 7-ray production region. 
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region filled with 
hot photon gas 




e + e ' 



Outside the blob - 
not boosted radiation 



Fig. 1. — A sketch of the model: a blob of proper radius -R blob (region filled with red color) 
moves with a bulk Lorentz factor T through a region of typical size R filled with a hot 
photon field (yellow region). Protons are accelerated and emit synchrotron radiation inside 
the blob with magnetic field strength B' . The produced synchrotron emission is assumed 
to be isotropic in the blob frame. The proton-synchrotron 7-rays can be absorbed due to 
pair production on the soft photon field. The pairs created in the blob produce detectable 
(Doppler boosted) synchrotron emission, while the emission of secondary pairs produced 
outside the blob is not Doppler boosted and therefore not detectable 
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2. Model description 

A sketch of the model adopted in this paper is shown in Fig. (TJ and the main ingredients 
of the model are described in the figure caption. 

2.1. Primary 7-rays 

The primary 7-rays are produced through synchrotron radiation of protons. Generally, 
in such a scenario the energy is stored in the magnetic field and episodically can be 



transferred to protons o 
particles (see for details 



extremely 



Aharonian 



high energies forming a non-thermal population of 



20001 ) . The energy released in non-thermal protons can 



be expressed through the strength of the magnetic field B' and the source radius -R(, lob : 

where k is the fraction of the blob magnetic energy transferred to accelerated protons (all 
physical quantities in the blob rest frame are primed). 

The proton synchrotron model for blazars works in the case of extremely effective 
acceleration of protons with an energy spectrum which continues up to ultra high energies, 
being limited either by the confinement in the accelerator (so-called "Hillas criterion"): 

S-SS-lO'^V^Uv, (2) 



or by synchrotron losses: 



B' 



-1/2 



E m <2-lV[ — ) TeV. (3) 



The synchrotron cooling time of protons 



E' \ - 1 / B' 
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is comparable to the proton confinement time assuming Bohm-type diffusion 

con " k d \WTeVj VlOOGy 1 Vl0 15 cmy' ' {) 
where kd is the ratio of the proton diffusion coefficient to the Bohm one. We note that the 
confinement time t' con cannot be shorter than light crossing time 

Css * 3 • 10 4 ( s. (6) 



10 15 cm / 

Given the identical dependencies of times described by Eqs.(jlI)-(jSJ) on proton energy, the 
cooling regime is defined by the following parameter: 

trnn 6 ( B \ ( R^cfo 



^ t^ yn k d V100G J Vl0 15 cmj ' ^ 
implying fast cooling for £ > 1 and slow cooling for £ < 1. In the case of slow cooling only 
a fraction £ of the proton non-thermal energy will be emitted through the synchrotron 
channel. 

In the blob frame the synchrotron emission is expected to be isotropic with spectral 
energy distribution (SED) extending up to the energy 

£maxH^0.2£TeV (8) 

if the proton maximum energy is given by Eq.Q the Hillas criterion Eq. ([2|) , i.e. £ < 3. 
In the case of large production region (i.e. £ 3> 1), the location of the SED maximum is 
determined by synchrotron losses and is expected to occur at 

E^-OATeV. (9) 

In the laboratory frame the proton emission of such a blob is characterized by 
luminosity of 

^ ^ (wd (lM S (i^g) 4 (lo4v) (s) 4 «*'- ■ < 10 > 
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where 5 is Doppler boosting factor, with typical variability time-scale of 

t var ~2-10 3 min(l,O (^-] (—^-z) (-^) s. (11) 



30 J Vl0 7 TeV/ ViooG^ 
Formally, the VHE spectrum of the boosted proton synchrotron may extend up to 

E max ~ 10 min(l, e/3) ( ^ ) TeV . (12) 



.30, 

However, we note that the actual shape of the spectrum close to the cutoff may be rather 
smooth, with a significant fraction of particles above the formal cutoff energy (e.g. for 



Zirakashvili fc Aharonian 



2003). This effect 



non-relativistic diffusive shock acceleration; see 
may significantly relax the constraints imposed by Eq. fll2jl . given the quadratic dependence 
of the synchrotron photon energy on the energy of the parent particle. 

For model calculations, in this paper we assume the spectrum of non-thermal protons 
to be a power law with exponential cutoff (N(E) = Nq E~ p exp(—E/ E c )), with the cutoff 
energy E c defined according to Eq.Q, i.e. we assume a very high acceleration efficiency. 
We consider two cases for the proton power-law index p: (i) the conventional value close to 
2, and(ii) a very hard case with index p = —0.5, as predicted by the converter mechanism 



f lDerishev et al 



20031 ) . The emission is assumed to be produced in the slow cooling regime. 
This approximation is valid for time intervals shorter than the variability time scale defined 
by Eq.f llip . Under this assumption, the VHE emission component is characterized by 
a photon index 1.5 in the case (i); and by the hardest possible photon index for the 
synchrotron radiation -namely 2/3- in the case of the converter mechanism (ii). We note 
that the latter case involves VHE spectra harder than conventionally accepted, but this 
case can be realized only in the slow cooling regime. Otherwise (i.e. in the fast cooling 
regime), the cooling mechanism should modify the proton distribution resulting in VHE 
spectra with photon index close to 1.5 for dominant synchrotron cooling; or close to 1 for 
dominant adiabatic losses. 
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2.2. Internal absorption 

There are several possible sources of UV and soft X-ray emission close to the base 
of the jet relate e.g. to accretion disk or corona. This emission may be reprocessed by 
matter surrounding the jet. This leads to the formation of the so-called Broad Line Regions 



(BLRs), which are characterized, in the case of powerful b 



luminosities of 10 45 erg s 1 (see e.g. 



Tavecchio &: Ghisellini 



azars, by a size of ~ 10 cm and 



20081 ) . Such dense photon fields 



imply significant 77 absorption, at least if the production region is located close to the jet 
base. The 77 optical depth is estimated as: 

r(£ 7 ) ~ 0.2a T Rn ph (3.5m 2 c 4 /£ 7 ) , (13) 

where <tt is the Thomson cross-section, R is the 7-ray travel distance in the photon field 
and n p h is the density of target photons. Since the size of the region filled by target-photons 
is larger than the travel distance R, the lower limit on the luminosity of the photon field, 
for the given optical depth r, is estimated as: 

L > 47rj R 2 ^P^ ~ 10 42 r ( ^ \ 1 ( -^—) ergs" 1 , (14) 

p h ~ 4 ViooGevy V 10 cm / 

where r is the maximum opacity which occurs for the 7-ray of energy E 1 . In general, the 
photon field required for the internal absorption scenario has a low luminosity and may be 
undetectable (it is not Doppler boosted!). In order to get an arbitrary hard spectrum after 
internal absorption, we assume the target photon field to be a gray body, i.e. a diluted 
Planckian distribution, characterized by the temperature T and the dilution coefficient (. 



2.3. Secondary emission 

The energy of the absorbed 7-ray is transferred to an electron-positron pair. Since 
the internal absorption scenario requires large optical depths r > 1, a significant fraction 
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of energy given by Eq.fJTJ goes to secondary electrons. The observational appearance of 
these secondaries depends strongly on the site of their production. Namely, if the pair is 
created outside the blob, emission of these electrons will not be boosted and thus remain 
undetectable. On the other hand, if the electrons are produced in the blob, they will be 
isotropized and emit synchrotron radiation due to the strong magnetic field in the blob. 
This radiation component can be detected because of strong Doppler boosting. 

In the blob reference frame, the target photon field is strongly anisotropic. Thus, the 
optical depth in the blob depends on the direction of the 7-ray with respect to bulk velocity. 
Since we assume the blob to be homogeneous, we introduce the optical depth r in , averaged 
over the 7-ray , to characterize the absorption in the 

blob. The corresponding values are shown in Table [TJ 

It is possible to estimate the optical depth r in in the blob since basically all the emission 
is focused towards the direction of the proper motion. Indeed, the optical depth for a 7-ray 
propagating in the direction of the proper motion will be 

nn(Ej) ~ 0.2a T i^ lob rn ph (3.5m 2 c 4 /£ 7 ) , (15) 



where T is blob bulk Lorentz factor (for detail see e.g. iBegelman et al.ll2008l ). Thus, a 
simple relation^ between size of the BLR region, size of the blob, maximum optical depth 



r, blob optical depth rj n and bulk Lorentz factor can be written as 

(16) 



r R 

This ratio indicates that in the case of a compact region filled with photon gas, R ~ 10 17 cm 



2 In our calculations the averaged optical depth is defined as e Tin =< e T >. 
3 This relation was used to derive the sizes of the blob listed in Table [Q 
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the optical depth in the blob is quite high: 

(s) Go^r • « i7) 

The injection spectrum of secondary electrons depends on the photon index of primary 
7-rays, target photon field, bulk Lorentz factor of the blob and the internal optical depth. 
If the target photon field is characterized by a peak energy e, then the maximum injection 
rate in the blob occurs at energy 

E' e ~ 5 ( ^— Y 1 ( — ^ GeV. (18) 
VIOeV/ V 20 / 

However, we have to note that depending on the slope of the primary 7-ray spectrum, this 
value can change significantly Since the synchrotron cooling time of these electrons, 

is very short (compared to both the typical time scales for the system and the Compton 
cooling time of electrons), 

4 " 7 ' 103 (iSv) " © _1 (to) ( 10W) (t) " 8 ' W 

the entire absorbed energy will be immediately released by secondary electrons through the 
synchrotron channel. 

In the case of large internal absorption or high bulk Lorentz factor, the secondary 
synchrotron component has a broad distribution centered at 

r\~ 2 /<5\ / e \- 2 / B' 

1 keV. (21) 



20y V30/ VIOeV/ ViOOG 

The variability time-scale of the synchrotron radiation of secondary pairs is determined 
by the change of the injection, i.e. by the change of primary 7-ray component. In the case 
of small internal opacity and assuming that protons are distributed over the energy interval 
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between 1 GeV and 10 7 TeV with E 2 -type spectrum, the luminosity of the secondary 
synchrotron radiation is estimated as 

L sec * 10- J^-) 3 ( *) V A V erg s- 1 . (22) 

In the framework of the discussed model an important relation can be established 
between the slopes of the intrinsic 7-ray spectrum and the highest energy part of the 
secondary synchrotron components. This part of the synchrotron spectrum is produced 
by electron-positron pairs which are created significantly above the threshold of the 77 
interaction, thus it is possible to use the asymptotic limit of the cross section. Since one of 
the secondary electrons receives almost the all parent 7-ray energy, the cross section can be 
approximated as 

Then, the spectrum of the secondary pairs, which is determined by the intrinsic spectrum 
of VHE 7-rays, is: 

dE e dt J 7 dE e dE y E e dE e ■ 1 } 

In particular, if the intrinsic 7-ray spectrum is a power law in this energy band, with a 

photon index s, then, since the dominant cooling mechanism is synchrotron radiation, the 

energy distribution of the secondary leptons is a power law with the index s + 2 and the 

high energy part of the synchrotron spectrum is described by a power law with photon 

index (s + 3)/2. We note that even for a very hard intrinsic 7-ray spectrum of s ~ 1.5, the 

synchrotron emission of secondary pairs will be characterized by a photon index ~ 2.25. 

Such behavior is expected at energies 

r \ ~V 5 \ ( e \ -2 / B' 



^ 200 {20) UJWJ UoGj UeV - (25 > 

For numerical calculations, we assumed the blob to be homogeneous. The pair 
production kernel, i.e. the energy distribution of secondary electrons produced by a 7-ray 
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of a certain energy, was calculated using anisotropic differential pair production cross 
section convolved with boosted Planckian distribution and averaged over the initial 7-ray 
direction. The injection rate of electrons was calculated by convolving the pair production 
kernel with proton synchrotron spectrum multiplied by the factor (1 — exp(— r- m )). The 
energy distribution of electrons was calculated using the approximation of continuous losses 
accounting for dominant synchrotron losses only. The secondary synchrotron emission was 
calculated using the obtained distribution of electrons. The variability properties of this 
radiation component are related to the variability of the intrinsic 7-rays as well as to the to 
the change of their absorption rate. 



3. Broadband spectra of 1ES 0229+200 and 3C 66A 



To demonstrate the potential of the proposed model for the explanation of very 
hard intrinsic 7-ray spectra, we focused on two distant objects, namely 1E S 0229+200 
(z = 0.1396) and 3C 66A [z = 0.444, though this value is debated, see e.g. 



Abdo et al. 



201ll ). These two BLLacs have different 7-ray properties. In particular, 1ES 0229+200 



shows VHE 7-ray emission without si gnificant flux or spectra l changes between two HESS 



measurements separated by one year ( lAharonian et al. 



20071 ). Moreover, Fermi LAT was 



not able to detect GeV emission from the direction of 1ES 0229+200. In the SED plot, 
the upper limit of the GeV flux appears below the TeV flux corrected for intergalactic 
absorption. 



The blazar 3C 66 A shows a variable VHE signal, as seen with VERITAS ( lAcciari et al. 



2009 



Abdo et al 



20111 ). with a 6% Crab flux flaring episode. Fermi LAT collaboration 



reported a significant GeV 7-ray excess from the sourc e. Moreover, an i ncrease of the GeV 



flux simultaneously with the VHE flare was observed ( lAbdo et al 



201ll ). Importantly, the 



GeV flux level exceeds significantly the de-absorbed VHE flux, thus a smooth connection 
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of these two radiation components with a single emission seems difficult to achieve, for the 
assumed redshift z=0.444. 

To study the impact of EBL on the VHE spectra, we corrected the reported 
7-ray spectra for i nterga lactic absorption using two versions of the EBL model by 



Franceschini et al. 



( 120081 1: (i) as in the original paper (FLO) and (ii) scaled up by a factor 



of 1.6 (F1.6). The 



Levenson fc Wright 



a fter c ase was considered in order to satisfy the lower limits claimed by 



( 120081 ). This simple treatment of the EBL and the related calculations 
of intergalactic absorption allows us to ignore many details of different EBL models, and 
focus on the main objective of this paper, namely the explanation of hard intrinsic 7-ray 
spectra in blazars. Note that the two EBL templates used here cover a broad range 
of different realizations of the EBL described by recent theoretical or phenomenological 
models, at least as long as it concerns the calculated optical depths. 

The optical depth for a high energy photon E 1 traveling through the intergalactic 
medium from a source at redshift z to the observer, taking into account the cosmological 
distance and the EBL evolution, is 

z 2 00 

T 11 (E y ,z) = c Jdz'— J dx-^ J den 1 (z',e)a yy (E 1 (l + z'),e,x) , (26) 

2m| e 4 

where ^7 is the cosmological line element; x — 1 — cos# is the angle between the interacting 
photons; n 1 is the number density of the EBL as a function of redshift and soft-photon 
energy; and <r 77 is the pair production cross section. In Fig. [2] the VHE 7-ray optical depths 
(left panel) and attenuation factors (right panel) for the two blazars are shown, for the two 
EBL levels: F1.0 (solid lines) and F1.6 (dashed lines). The calculated attenuation was used 
to reconstruct the initi al spectra from the observed da ta by H.E.S . S. on 1ES 0229+200 



( lAharonian et al 



20071 ) and by VERITAS on 3C 66A (lAbdo et al 



20111 ). The resulting 



spectra are shown in Fig. [3] for 1ES 0229+200, and in Fig. H]for 3C 66A. In both figures 
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black points correspond to observed data, red points to the spectra reconstructed with the 
F1.0 EBL model, and blue points to the spectra reconstructed with the F1.6 EBL model. 

The reconstructed spectra are significantly harder compared to the observed ones. 
In particular, in case of high EBL flux (F1.6), the spectra have r; nt ~ 1 and 0.2 for 
1ES 0229+200 and 3C 66A, respectively; i.e. they would be significantly harder than the 
conventional value of 1.5. 



3.1. The case of 1ES 0229+200 



To study the case of 1ES 0229+200, we have combined the reconstr ucted VHE data 



with archive X-ray and 



( jCostamante et al. 



optical data from SWIFT ( iTavecchio et a. 



2009) and 



BeppoSAX 



20091 ). The 



20021 ). together with Fermi LAT observations f lAbdo et all 
observational data are summarized in Fig. [3j We have applied the internal absorption 
scenario as described in Section [2] to reproduce the VHE spectrum together with X-ray 
spectrum, for both levels of intergalactic absorption and considering both indices of the 
power-law proton distribution (p ~ 2 and p = —0.5). 

In the case of a soft energy distribution of protons (p > 2), the flux upper limit 
obtained with Fermi requires a lower-energy cutoff (E[ c ) in the proton energy spectrum at 
very high energies, with both levels of EBL absorption. The exact location of this cutoff 
depends on different model parameters, e.g. intrinsic optical depth, Doppler boosting factor 
etc., but to satisfy the Fermi upper limits the resulting 7-ray spectral break should occur 
close to ~ 100 GeV. Thus, the cutoff in the proton spectrum should be located roughly at 



E> ~ 2 ■ 10 6 ( 
lc 1 100 G 



-1/2 



5_ 
30 



-1/2 



TeV 



(27) 



This value is very close to the highest possible energy of the accelerated protons, thus 
the Fermi upper limits basically exclude the possibility of a proton spectrum significantly 
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steeper than dN/dE oc E 2 . 

Since the VHE spectrum obtained from 1ES 0229+200 seems to show no significant 
changes on a yearly time scale, i.e. on a time scale much longer than the one defined by the 
cooling time (see Eq. U}, the proton spectrum is expected to be steady. A steady proton 
distribution with power-law index p = 2 can be formed in two different ways: (i) with an 
almost mono-energetic continuous proton injection (e.g. through converter mechanism) in 
the fast synchrotron cooling regime; and (ii) with a conventional acceleration spectrum in 
the slow cooling regime. A very hard steady proton distribution with p = —0.5 requires can 
be formed in the slow cooling regime when an acceleration mechanism similar to converter 
mechanism is responsible for the particle acceleration. 

The resulting model parameters are summarized in Table [T] and the corresponding 
curves (Fits 1-3) are shown in Fig. [3j 

3.1.1. F1.0EBL level 

In the case of the EBL level F1.0, the de-absorbed VHE spectrum has a photon index 
close to r int ~ 1.5. For a proton distribution with p = 2, the proton synchrotron radiation 
below the peak has a photon index close to 1.5, i.e. formally it can explain the VHE 
data points without invoking internal absorption. Thus, in this case the key question is 
whether the internal absorption scenario can provide a consistent explanation of the X-ray 
component. 

Given the strict upper limits provided by Fermi, which are at the level of the 
extrapolation in the HE band of the r int ~ 1.5 VHE spectrum, the available energy 
budget for the secondary pairs is quite limited, unless a higher emission can be effectively 
suppressed in the Fermi-LAT band. This could be achieved either by assuming a broad 
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energy distribution of target photons extending to X-ray energies, so to provide a significant 
attenuation also in the GeV band, or by introducing a very high lower-energy cutoff in the 
proton distribution. 

In absence of these two conditions, the X-ray synchrotron flux of the secondary pairs 
would be approximately an order of magnitude below the reported X-ray fluxes. Therefore, 
in this specific case, the internal absorption scenario requires additional ad-hoc assumptions 
to provide a self-consistent interpretation of the TeV and X-ray data. 

These additional assumptions instead are not needed in the case of a hard proton 
spectrum (p = —0.5). The latter can provide both the energy budget to explain the 
X-ray data and GeV fluxes below the Fermi-LAT limits, as shown in Fig. [3] (Fit 2, whose 
corresponding parameters are given in Table [T]) . 



3.1.2. F1.6EBL level 

In the case of high EBL flux (model F1.6), the de-absorbed VHE spectrum has a 
photon index close to r int ~ 1, i.e. harder than the unabsorbed synchrotron spectrum from 
a proton distribution with index p ~ 2. Internal absorption allows the hardening of the 
TeV spectrum to the required level, but in the case of a conventional proton distribution 
the discrepancy with the Fermi-LAT upper limits is very strong. To avoid the conflict with 
Fermi-LAT data we need again to suppress the GeV emission, by introducing additional 
assumptions such as an effective absorption of GeV 7-rays (e.g. by X-rays) or a very high 
lower-energy cutoff (at 10 6 TeV in proton energy). However, these assumptions can hardly 
be endorsed without an additional observational or theoretical justification. 



On the other hand, a very 



mechanism (IDerishev et al. 



2003|) 



rard proton distribution as predicted by the converter 
can accommodate the Fermi LAT upper limits. Even so, 
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the synchrotron spectrum would not be sufficiently hard to explain the 7-ray spectrum 
corrected for the higher EBL flux (model F1.6). In this case the internal absorption becomes 
a key requirement to further harden the initial proton synchrotron spectrum (see Fits 1 and 
3 in Fig. [3j and the corresponding parameters in Tabled]). 

The synchrotron radiation of secondary electron-positron pairs, calculated self- 
consistently with the hard VHE component, can explain the X-ray flux of 1ES 0229+200, 
with the caveat that the X-ray data are not simultaneous with the TeV observations. A 
characteristic feature of the secondary synchrotron radiation is its broad spectral extension 
up to hard X-rays. This prediction can be tested with hard X-ray instruments like Suzaku, 
or with future missions NuStar and Astro-H. 



3.2. The case of 3C 66A 



For the blazar 3C 66A, we have combined the VHE spectrum reported by VERITAS 
with the spectrum detected by Fermi-LAT during the VHE fl are, together with the available 



X-ray/optical data from MDM and Swift (lAbdo et al. 



20111 ). The observational data are 



summarized in Fig. H] and EJ We have applied the internal absorption scenario as described 
in Section [2] to fit the VHE spectrum, again considering both levels of EBL absorption 
(F1.0 and F1.6). 

To correct for intergalactic absorption, the redshift for the source we adopted is the 
one most often cited and used in the literature, z = 0.444. It should be noted, however, 
that the redshift of this source is not yet firml y established, and thus one cannot exclude 



Prandini et al 



( 120 id ) suggested that the 



that the source is located closer. In particular, 
redshift should not exceed 0.34. This conclusion is based on the belief that the initial VHE 
7-ray spectrum cannot be harder than the GeV spectrum measured with Fermi LAT. On 
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the other hand, if the redshift is indeed > 0.4, the TeV and GeV parts look quite different, 
and not part of a single component. Even though, this does not imply that they must be 
of different origin. In fact, our model can explain both components with a single proton 
population, as parts of the smooth proton synchrotron spectrum which is then deformed by 
the energy-dependent internal absorption. 

3.2.1. F1.0EBL level 

In this case the de-absorbed TeV spectrum is rather flat, with photon index r int ~ 1.7, 
while the HE component is characterized by a similar photon index T ~ 1.8 but with at 
higher flux. A good agreement between the GeV and TeV spectra can be achieved assuming 
a proton energy distribution with power-law index p = 2. 

A weak internal absorption (with maximum optical depth of about r = 1.6) allows 
modification of the VHE spectrum to the required photon index (Fit 5 in Fig. H]), while the 
HE part is reproduced by the unmodified synchrotron spectrum. The synchrotron emission 
of secondary pairs can explain the X-ray spectrum obtained with Swift but not the optical 
MDM data, which require an additional radiation component. 

The physical parameters used in this model may appear quite extreme (see Table HJ 
Fit 5). In particular, the very small value of the Doppler factor has been chosen to avoid 
7-ray excess above 1 TeV, and this consequently leads to a dramatic increase of the required 
energy budget. In fact, there is a more natural way to suppress the flux level above 1 TeV, 
namely assuming a less efficient acceleration process. In this way the Doppler factor and 
B-field may be increased, while the required energy budget will be significantly reduced. 
The detail study of this possibility will be discussed elsewhere. 
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3.2.2. F1.6EBL level 

For the high EBL flux a very small photon index of r; nt ~ 0.2 is required. Remarkably, 
even such an unusual photon spectrum can be explained by internal absorption with a 
higher target photon temperature and slightly larger optical depth (see Fig. HI Fit 4). With 
a certain combination of model parameters, the flux of the synchrotron radiation from 
secondary electrons can match the levels detected in the optical band, as is demonstrated 
in Fig. \5\ Fits 4a and 4b. For a small production region, the main fraction of the secondary 
pairs are produced outside the blob. Their radiation is not Doppler boosted and, therefore 
cannot be detected. For the pairs produced inside the blob, the secondary synchrotron 
radiation is Doppler boosted and thus it can contribute significantly to the observed 
fluxes. We note, however, that for this source we did not succeed to find a combination of 
parameters which could explain both the optical and X-ray fluxes together, by synchrotron 
radiation of secondary electrons. Since the internal absorption scenario requires a significant 
attenuation of the VHE radiation over approximately two decades (see Figs. [3] and EJ), 
the secondary synchrotron component has to be at least 4 decades broad (with additional 
broadening related to the relativistic motion of the production region). The strong magnetic 
field required in the proton synchrotron model provides fast cooling of the pairs, thus the 
radiation spectrum will be featureless, without a cooling break. In the case of a small radius 
of the production region, the effective particle injection in the blob may be rather narrow. 
But in this case the flux level would be significantly below the observed flux. 
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Table 1. The combination of the parameters used for the calculations of the model curves 
in Figs. |3l H] and [5] for the different EBL levels (first row): p is the power-law index of the 
proton distribution, B' is the magnetic field inside the blob, T is the temperature of the soft 
photon field, r is the optical depth for the entire source of soft photons in a region of radius 
R, ri n is the optical depth inside the blob, -Rb lob is the proper radius of the blob, T the bulk 
Lorentz factor, 5 the Doppler factor, L p h the luminosity of the soft photon source, L'^ is the 
intrinsic luminosity of the 7-ray source before Doppler boosting and internal absorption. 



1ES 0229+200 3C 66A 
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Fit 2 


Fit 3 


Fit 4a 


Fit 4b 


Fit 5 


EBL 


F1.6 


F1.0 


F1.6 


F1.6 


F1.6 


F1.0 


V 


-0.5 


-0.5 


-0.5 


2.2 


2.2 


2 


B'(G) 


80 


40 


100 


100 


100 


1.2 


T(K) 


7 x 10 3 


5 x 10 3 


10 5 


8 x 10 4 


8 x 10 4 


5 x 10 4 


T 


3 


3 


5 


2 


2 


1.6 


Tin 


0.9 


0.8 


1.2 


0.6 


6 x 10~ 3 


0.26 


R'biob (cm) 


10 15 


5 x 10 15 


5 x 10 15 


5 x 10 17 


5 x 10 15 


10 18 


R (cm) 


3 x 10 16 


6 x 10 17 


2.1 x 10 17 


7 x 10 19 


7 x 10 19 


6 x 10 19 


r 


10 


30 


10 


40 


40 


10 


5 


11 


8 


8 


15 


15 


4 


L ph (ergs -1 ) 


2.3 x 10 41 


2.8 x 10 42 


3 x 10 43 


4 x 10 45 


4 x 10 45 


2 x 10 45 


L'^ (ergs -1 ) 


5 x 10 40 


1.6 x 10 41 


2.9 x 10 41 


7 x 10 42 


7 x 10 42 


4 x 10 44 
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Fig. 2. — Left panel: The optical depth r for 7 rays interacting with the EBL. The upper 
set of lines is for z = 0.1396 (distance of 1ES 0229+200) while the lower set is for z = 0.444 
(distance of 3C 66 A). The dashed lines are calculated for the EBL model F1.6, the solid 
lines for the EBL model F1.0. Right panel: The attenuation factors e~ T corresponding to 
the optical depths shown in the left panel. 
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Fig. 3. — The overall SED of the blazar 1ES 0229+200, together with model curves calcu- 
lated with the proton synchrotron scenario plus internal absorption. The set of parameters 
used in the calculation of Fits 1 (green), 2 (red) and 3 (blue) are presented in Table [D Inner 



panel: the VHE spectrum as observed by H.E.S.S. (lAharonian et al. 



20071 . black points). 



The same spectrum corrected for EBL absorption with model F1.0 (red points) results in an 
intrinsic power-law index of r int ~ 1.5, while using the higher EBL model F1.6 (blue points) 
it yields T int ~ 1.0. Outer panel: The same data as in the inner panel with the addition 
of the SWIFT data (X-ray and optical bands, magenta points) and of the BeppoSAX data 
(X-ray band, blue points). The solid black lin es between 30 Me V and 30 GeV are the up- 



per limits based on FERMI LAT observations ( lAbdo et al. 



20091 ). and calculated assuming 



power-law 7-ray spectra with photon indices r int = 2 (a), r int = 2.5 (b) and r int = 1 (c) 
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Fig. 4.— The VERITAS spectrum of 3C 66A (jAbdo et al. 



2011 



flare dataset, black points). 



The same data corrected for intergalactic absorption using the EBL model F1.0 (red points) 



result in a power-law intrinsic spectrum (dN/dE oc E~ r ' int ) with photon index Tint ~ 1.7. 
Using the EBL model F1.6 (blue points), instead, the photon index becomes r int ~ 0.2. 
The two model lines (labeled 4 and 5) are calculated using the parameters presented in 
Tabled] (Fit 4 and 5). Inner panel: Zoom out of the plot to include the GeV band. The 



data points 



VERITAS flAbdo et al. 



corresponc 



to t he FERMI LAT observations performed simultaneously with 



20111 ). 
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Fig. 5.— The overall SED of 3C 66A as observed in VHE (VERITAS: 



Abdo et a 



GeV (FERMI LAT) 7-rays together with optical (MDM) and X-ray (Swift) data flAbdo et al. 



201 lh and 



201ll ). The curves represent the model calculations performed for the combinations of pa- 
rameters reported in Table [TJ Fit 4a and 4b differ only in the size of the 7-ray production 
region (large or small, respectively), i.e. if the majority of the electron-positron pairs are 
produced inside or outside the relativistically moving blob. The dashed lines 4 and 5 show 
the spectra before internal absorption. 
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4. Summary 



One of the most challenging issues of the physics of TeV blazars is the rather hard 
intrinsic 7-ray spectra of some representatives of this high-energy source population. 
Actually, the reported spectra themselves are steep, with photon indexes r int > 3. However, 
the spectra of VHE 7-rays after being corrected for the energy- dependent intergalactic 
abso rption become very hard, in some cases as hard as r; nt = 1.5, assuming the EBL flux 



from 



Franceschini et al. 



Dominguez et al. 



( 120081 ). Note however that some other recent models, in particular 
( 120101 ). give very similar results. Explanation of such spectra faces serious 
difficulties within the standard blazar models. Moreover, even a slight increase of the EBL 
flux at optical and near IR wavelengths compared to the benchmark models (which, given 
the significant uncertainties in the derivation of the EBL fluxes, cannot be excluded) results 
in unusually hard intrinsic spectra with r int < 1. 

In this paper we studied the applicability of the idea of internal absorption of 7-rays 
produced by highly magnetized blobs as a result of synchrotron radiation of protons. While 
the main aspects of the model have been developed and discussed in our previous paper 



( lAharonian et al 



20081 ). in this work we tried to understand whether the model can be 
applied to specific objects. For this purpose, we have chosen two "difficult" representatives 
of this source population, 1ES 0229+200 and 3C 66A, and allowed the EBL flux to be 



somewhat higher than the lowest possible fluxe s. In particular, for t 



le EBL level consistent 



Levenson &: Wright 



(120081 ) . the intergalactic 



with the lower limit around 3.6 /im claimed by 
7-7 de-absorption of the VHE flux detected from 3C 66A during a flaring episode results 
in an extremely hard spectrum, with photon index Tint = 0.2. While such a hard spectrum 
cannot fit into any existing VHE 7-ray production model, the scenario of internal absorption 
of 7-rays produced via synchrotron radiation of protons provides a reasonable explanation 
of both the hard TeV spectrum and the high GeV flux observed during the flare. 



-30 - 



In the case of 1ES 0229+200, the internal absorption scenario calls for an extremely 
hard proton distribut ion, which can be provided for example by the "converter mechanism" 



f lDerishev et al 



20031 ). The synchrotron radiation of secondary electron-positron pairs 



produced inside the blob results in an additional (Doppler boosted) radiation component, 
which can provide a self-consistent interpretation of the non-thermal X-ray emission in both 
objects. However, since in the framework of the internal absorption scenario the secondary 
synchrotron component is expected to be quite broad and featureless, we failed to find a 
set of parameters which explains simultaneously both the optical and X-ray data obtained 
from 3C66A. 

The intrinsic absorption scenario allows a natural explanation of the very hard 
intrinsic TeV spectra at the cost of a large attenuation of the 7-ray flux around 100 GeV. 
Quite remarkably, despite the significant attenuation, this scenario does not enhance too 
dramatically the required energy budget. Indeed, since the correction for intergalactic 
absorption requires the initial 7-ray spectra to be very hard to begin with, with a photon 
index < 2, the energy requirement to reproduce such a spectrum is determined by the 
highest energy part of the spectrum, which is not affected by internal absorption. Typically, 
the enhancement of the energy budget introduced by internal absorption does not exceed 
a factor of 5, which can be easily compensated by a slightly enhanced Doppler boosting of 
the radiation. Therefore the required intrinsic 7-ray luminosities remain quite modest (see 
Table [TJ). Given the very high (almost 100%) efficiency of conversion of the proton energy 
to 7-rays through the synchrotron radiation in the magnetized blob, the suggested scenario 
can be treated as quite effective. 

It is apparent that the interpretation of very hard 7-ray spectra by the internal 
absorption depends on the "right" choice of several model parameters, especially if we want 
to explain the X-ray data by synchrotron radiation of secondary electrons. Nevertheless, 
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we should note that the level of hardening of the VHE 7-ray spectrum depends, in fact, 
only on the temperature of the photon field and the optical depth. On the other hand, 
while the flux ratio of the VHE 7- and X-ray components depends mainly on the size of the 
production area, the location of the secondary synchrotron peak is sensitive to the strength 
of the magnetic field, to the photon temperature and the bulk Lorentz factor. Finally, we 
note that certain radiation features of the scenario do not depend on the model parameters 
at all. This concerns, in particular, the slope of the 7-ray spectrum at GeV energies, and 
the shape of the X-ray spectrum. 

The dependence of the results of radiation properties on several parameters limits, 
to a certain extent, the predictive power of the suggested model. This is a consequence 
of complex environment in blazars where several radiation and absorption processes can 
proceed simultaneously. In this regard, the often used one-zone models with consideration 
of only synchrotron and IC radiation components produced in the same region, are quite 
useful for understanding the basic aspects of the problem, but can hardly properly describe 
the complex scenarios that take place in blazars. In particular, the results of this paper 
demonstrate that the internal absorption not only cannot be a priory excluded from 
the consideration, but, in fact, in some cases can be invoked for better explanation of 
observations of TeV 7-ray blazars. 

The authors would like to thank Dr. L. C. Reyes for making available the VERITAS 
data of 3C 66A flare and Prof. F. Tavecchio for the SWIFT data of 1ES 0229+200. 
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